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We derive the equation of matter density perturbations on sub-horizon scales around a 
flat Friedmann-Lemaitre-Robertson- Walker background for the general Lagrangian density 
fiR, Q) that is a function of a Ricci scalar R and a Gauss-Bonnet term Q. We find that the 
effective gravitational constant generically scales as distance squared at small distances. The 
effect of this diminishing of the gravitational constant might be important in the gravitational 
dynamics of cosmic objects such as galaxies, which can be in principle tested by observations. 
We also provide the general expressions for the effective anisotropic stress, which is useful 
to constrain modified gravity models from observations of large-scale structure and weak 
lensing. We also find that there is a special class of theories which evade this unusual 
behaviour and that the condition to belong to this special class is exactly the same as 
the one for not having super-luminal modes with propagation speed proportional to their 
wavenumber. 

§1. Introduction 

Modifying the law of gravity is a well-known possible explanation for the origin 
of the accelerated expansion of the universe (for a general review of the acceleration 
mechanism, see e.g^^). Precise observations of cosmic structures are now becom- 
ing powerful methods to distinguish modified gravity scenarios from other scenarios 
such as a cosmological constant. The basic observation which enables us to hope it is 
indeed feasible to recognize among these theories the right one, is that modified grav- 
ity theories generically predict characteristic scale-dependent growth of the matter 
density perturbations, which thus leaves distinct patterns of the cosmic structures. 
In other words, the dynamics of matter perturbation becomes scale dependent, and 
each scale feels in general a different effective Newtonian constant, i.e. a different 
gravity law. Therefore by precisely measuring the patterns of the cosmic structures, 
we can in principle accept or rule out the modified gravity scenarios as the origin 
of the accelerated expansion. There have been a number of recent studies about 
the evolution of density perturbations for various types of the modified gravity the- 
ories, e.g. for f{R) and scalar-tensor theorie^^^IZli for R -\- f{Q) theories for 
fiR^Ry^uR^"" ,UR) theories,l23 and for the DGP-inspired Galileon field.EB 

In,12^ we have studied the basic properties of propagation of the perturbations 
around a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background for the 
general f{R,Q) theories with a perfect fluid. There, we found that there are four 
independent modes for the perturbations. Two of which have dispersion relations 
like ^2 = ASk'^ + Oik"^), where A is a background dependent quantity and ^ is a 
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determinant of the Hessian of f{R,G): 

" ~ dm dg^ \dRdg J ' ^ ^ 

evaluated for the background metric. For the special theories such as f{R) and 
R+f{Q) that give vanishing S, the leading term in the dispersion relation starts from 
order. On the other hand, the remaining two have the standard ones = c^k'^, 
where = P/p is the sound speed of the fluid. The former generally results in super- 
luminal propagation for short wavelength modes, while the latter is the standard 
dispersion relation for a perfect fluid. 

Based on the basic studies in,^^ in this paper we restrict ourselves to the non- 
relativistic matter (c<j = 0) and study in detail the properties of the matter density 
perturbation for f{R,Q) theories to derive the consequences that are linked with 
observations. We will provide general formulae for the parameters that are now 
widely used in literature to parametrize the deviation from GR for f{R, Q) theories. 
Interestingly, it turns out that the effective gravitational constant generically decays 
at short distance r, obeying a scaling oc r^. We also discuss a class of theories which 
evade this diminishing of the gravitational constant. 

§2. Modified Gravity Models 

We consider the f{R,Q) theories where cold dark matter(CDM) and baryons 
are minimally coupled to gravity. The corresponding action is given by 

S = j d^x^g fiR,g) + I d^x^g (2-1) 

where Gn is the Newton's constant, R is the Ricci scalar and Q is the so-called 
Gauss-Bonnet term defined by 

g = R^ - AR^.R^'^ + R^^^pR^^'^^ , (2-2) 

and Cm is the Lagrangian density for the CDM and baryons. 
The equations of motion for g^^^ are given by 

Here E^i, is the effective energy momentum tensor defined by 
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^^u = ;^(V^V,F - g^^DF + 2i?V^V,e - ^g^^RO^ - AR^^V^V,^ - AR.^V^V^^ 
+AR^,a^ + Ag^uR'^^VaVf,^ + 4i?^,/3.V"V'^C - i g^,V) , (2-4) 
where F, ^ and V are defined by 

V = RF + g^- f. (2-6) 
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T^y is the energy-momentum tensor for CDM and baryons. It is given by 



(2-7) 



where h/c represents baryons/CDM respectively. As usual, we impose a normaliza- 
tion condition u'ff'u^^^^ = — 1 for each component. 

2.1. Background dynamics 

We assume that the background spacetime is a flat Friedmann-Lemaitre- Robertson- 
Walker (FLRW) universe whose metric is given by 



ds^ = -dt^ + a^{t)6ijdx'dxK 
Then we find that the background equations are given by 
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-SirGNprn + HF- 2HF + 4H^^ - 8HH^ - 4H^^, 



(2-8) 

(2-9) 
(2-10) 



where pm — 

pib) + p{c) is 

the total matter density. 

2.2. Scalar-type Perturbation 

We consider linear scalar- type perturbations around the metric Eq. (|2-8p . We 
always work in the Newtonian gauge: 



-(1 + 2<F) dt^ + a^{t){l + 2^)5ij dx' dx^ . 



(2-11) 



We also decompose the energy-momentum tensor for the CDM and baryons as 

piA) ^ p{A) + ^p(A) ^ ^ ^(A)^^ ^(A) ^ g^yiA)^ ^ ^^^^^ (2-12) 

All perturbation variables for the matters are also defined in the Newtonian gauge. 
Then the perturbation equations, in Fourier space, are given by 



1 



SH^'F - -^<P - 3H<P = — 
a2 2F 



3(F + l2H'^C)d + S^H^^ - 6H{F + 8H'^(,)F 



+ 3H{6F + 4H'60 + ( -3(F^ + H) + ] {6F + 4H'6^) 



<P-HF = — 
2F 



F 

r2^'t I A u-is 



+ 



H6F -5F- m'^Si + AH-^i + (F + l2H'^i)F - 8Hi<P 
F ' 



^-\-F 
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AHi^ + 4^^ + 5F + A{H^ + H)dC 



(2-13) 

(2-14) 
(2-15) 
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5^"' + 3<? 



i-2 



0, 
0, 



(2-16) 

(2-17) 

(2-18) 
(2-19) 



The perturbation quantities 5F and 5^ are related to perturbations of Ricci scalar 
and the Gauss-Bonnet term through the equations: 



6F =Fr6R + Fg6g, 
6C =Fg6R + ^g6g, 



(2-20) 
(2-21) 



where = dF/dR and so on. In the second equation, we have used an identity 
= Fg. By perturbing the scalars R and Q, we have the following relations: 



5R 

sg 
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- 3^ - 12H^ + 3H'F 



(2-22) 



u2 

-2^{H^ + H)<P 



n2 



(2-23) 
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3.1. Catch-up of haryons 

Because the matters are minimally coupled to gravity, the baryon density per- 
turbation starts catching up CDM after the recombination. To explain this behavior 
in more detail, let us introduce a quantity A = 6^^^ — (5^^^ By using ()2T6p - ()2T9p . 
we can derive a following evolution equation for A, 



A + 2HA = 0. 

The general solution of this differential equation is given by 

dt' 



A{t) = a^(ti)Z\(ti 



(3-1) 



(3-2) 



where ti is an arbitrary constant and the two initial conditions A(ti) and A{ti) may 
be in general k dependent. If we take ti as a time just after the decoupling of CMB 
photons from baryons, the expansion of the universe can be well approximated by 
the one for the matter dominated universe, i.e., a{t) oc t^/^. Then the integral on 
the RHS of Eq. ()3-2p is mostly determined by its lower limit which becomes 



a'^{ti)A{ti] 



dt' 



O 



Ajti 
H{ti 



(3-3) 
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Since both matter and baryon density perturbations can grow only on the cosmo- 
logical time scale ~ 1/H{ti), the first term in (j3-2p remains of the same order of 
magnitude as the second one even at late times. Meanwhile the matter density per- 
turbations themselves grow from tiny fluctuations to more than to make cosmic 
structures. Therefore, a relative difference A/ 6^'^^ quickly decays and the baryon den- 
sity perturbation becomes almost the same as the CDM density perturbation. Note 
that we have used only the conservation equations for the matter fields to derive 
this result. The only place where the effect of modification of gravity enters is the 
evolution of the scale factor appearing in the integral in (|3-2p due to the modification 
of the Priedmann equation. However, since the effects of modification of gravity are 
negligible at high cosmological redshift, we expect that the evolution of A to change 
only slightly from the one in GR. 

At late times, when the modification of gravity becomes important, we have 
(5^'') ^ 6^'^^ . Then it is a good approximation to consider a total density perturbation 
6m defined by 

^^^^ 

Pm 

instead of considering d^'^^ and 6^^^ independently. The approximation 6m = 5^^^ = 
6^^^ yields another relation V^^^ = y*-^-* = Vm, which can be derived by substracting 
(|2T6D from (I2T71) . 

To conclude, baryon density perturbation catches up the CDM density per- 
turbation in the same way as in GR. Correspondingly, we can safely replace the 
conservation equations (|2-16|) - (|2-19|) by 

^2 

5m + S4>- -^Vm = 0, (3-5) 

Vm + ^ = 0, (3-6) 

to study structure formation at low redshift universe. In the following analysis, we 
use 6m and Vm for the matter density perturbation. 

3.2. Effective gravitational constant 

One generic feature of the modified gravity models is that the effective gravita- 
tional constant has a scale dependence over cosmological scales. In the literature it 
is common to define the effective gravitational constant by a Poisson equation 0: 

— ^ = -A7TGNQ{k,t)pm6m, (3-7) 

where we parametrize the deviation from GR by Q{k,t) (as in GR we have Q = 1). 
The modification of the Poisson equation will affect the evolution of the matter 
density perturbations, which can be tested by cosmological observations such as 



*■* Some literature (e.g.'^''^) define the eilective gravitational constant by using <1> instead of 
'I'. If we use <5, the effective gravitational constant then represents how strongly the space is curved 
by the presence of matter. To go from one to the another, we have to multiply or divide by a factor 



6 



Antonio De Felice and Teruaki Suyama 



weak lensing. We use ^, the time component of the metric perturbation, to define 
the effective gravitational constant. Therefore the gravitational constant obtained 
by measuring a force between two mass points is given by GnQ. 

In this subsection, we derive the general expression of Q{k, t) for f{R,Q) theories 
on sub-horizon scales which are relevant to cosmological observations. To this end, 
following the approach in,'3'Sl''121i:l2Q) ^jjj ^ subhorizon approximation under 
which we assume that time derivative of any perturbation variable is much less than 
its spatial derivative or equivalently, in terms of Fourier space, a term multiplied by 
k. Therefore if there are time derivatives of a perturbation variable in addition to 
terms of the same variable multiplied by k"^ in a perturbation equation, we drop the 
first from the equation. By this procedure, we are dropping the two fast oscillating 
modes among the four whose dispersion relation is given by cj^ = ASk^ {AS must 
be positive to avoid an instability of the FLRW universe) and are taking into account 
only the remaining two modes which evolve on cosmological time scale. 

With the sub-horizon approximation, the perturbation equations are given by 

-%f^J^{mi.p + SF + iHHi]-'-!^^^i!^ (3.8) 



a? 2Fa? \ ^ F 
AHi<F + + 6F + 4(i?2 + H)6(, 



'P + ^ = -j 



(3-9) 



i,2 

5R = 2^(2<P + 'F), (3-10) 

k^ " 
6G = S^ 



(3-11) 



J^2 

5m - -^Vm = 0, (3-12) 

Vm + 'F = 0, (3-13) 

where we again use the relations (j2-20p and ()2-2ip to relate {6F,6S,) with {6R,6Q). 
In Eq. (|3-8p we have neglected the contributions of the 'F field as it is of the same 
order of ^. This assumption turns out to be consistent when we write, as we shall 
see later on, ^ in terms of ^. We do not include ()2T4p as it is not necessary for our 
present purpose. 

To derive the effective gravitational constant, we first combine ()2-20p . (|2-2ip . 
(j3-10p and ()3Tip to express 6F and 5^ in terms of 'P and F. Substituting these 
results into (|3-9p . we obtain an equation which only contains <P and F. By solving 
the obtained equation with respect to <P, F is expressed solely by F, and at this point 
we can check the consistency of the assumption we considered earlier. By putting 
this relation back into equations for 5F and to eliminate F, both 6F and (5^ are 
also expressed solely by F. Then by combining these results with (|3-8p . we find that 
6m can be also written solely by F, from which we can derive the desired expression 
for Q{k,t): 

^ A^ + MmaHr „ 

By+B2lk/{aH)f + B;ilk/{aH)]*' 
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where time-dependent coefficients are given by 

Ai = F + 4|, (3-15) 



^9 = AH 



(3-16) 



Fr + %Fg{H' + H) + lQig{H' + H) 

Bi = {F + mif, (3-17) 
B2 = H\3F + 16Hi - 4i)[FR + 4Fg{2H^ + H) + 16H^(g{H^ + H)] 

+ 4HH\F + 4i){Fg + 4H^(g), (3-18) 

Bs = UH^H^{F^ - Fnig) = -UH^H^S . (3-19) 

In Ai and Bi, there are other terms that are functions of Fji,Fg and ^g. However, 
they are subdominant compared to F terms and we have neglected them in the 
above equations. We find that B^ is proportional to S. Therefore, it is absent 
for the special models where S = 0. As we will see, the behavior of the effective 
gravitational constant crucially depends on if the B^ is zero or not. 

Let us consider f{R) theories just to check that our result correctly reproduces 
Q given in the literature. In this case, C = ^ = -Fg = ^g = and we have 

F + AFnH^k/iaH)? 
F-^ + 3FFRH^[k/iaH)f 

This coincides with the one given in^^ 

3.3. Interpretation of the result 

Now let us consider (|3-14|) where we find some new interesting consequences. 
For long wavelength modes (A; — )• 0), we find that Q becomes independent of k: 

Q-# = ^, (3-21) 

where 

4^ = ^ + _ and Qtt = hF + '^Hi , (3-22) 

are respectively the squared speed and the coefficient of the kinetic term for the 
gravitational wave modes of the theory (se^^S}^^ Both these quantities must be 
positive in order to avoid ghosts and small scale Laplacian instabilities {c^rp < 0). A 
consequence of Eq. (|3-2ip is that an opposite sign for Qtt would lead to an effective 
repulsive gravity law. For short wavelength modes {k — )• oo), we find that Q is given 
by 

Q«^^. (3.23) 

We see that Q is inversely proportional to A:^, which means that gravitational 
constant decays on small scales. This is quite different from /(i?) theories be- 
cause in f{R), Q approaches a value 4/(3F) on small scales. Let us just refei!^^ 
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which found that /c~^ scahng of the effective gravitational constant is present for 
/(i?, R^uR^'^ ■, ^R) theories. This theory does not overlap with /(i?, Q) theories con- 
sidered in this paper as a quantity R^uapR^'^"^^ does not enter in their action. Fur- 
thermore, that theory in general may contain spin-2 ghost degrees of freedom. 

To see what is going on in real space, let us consider the gravitational potential 
sourced by a local distribution of matter density on the cosmological background. As 
the simplest case, let us consider a situation that a mass point is put at the origin, 
i.e., 5pm{x) = M5{x) 0- Then the solution of the equation (|3-7p in real space is 
given by 

mx t) = r A,+A,[k/{aH)f 

Jo k Bi+B2[k/[aH)f + B-i[k/[aH)t ' " 

where r = Although we can do the integration analytically, the result is rather 
complicated. Here we decide to provide only asymptotic behavior of (|3-24|) at large 
and small distances from the source. 

For large r, the dominant contribution to the integration comes from a region 
near A; = 0. Then we can drop the higher order terms in k appearing in the integrand 
and we have 



l^(x,t) = -^^-^ + 4 . (3-25) 



GnM _4j_ 1 
■ 2Qtt 

As expected, the gravitational potential is inversely proportional to r with its am- 
plitude multiplied by background quantities. On the other hand, for short distance, 
the higher order terms in k become important and we have 

= const. + r + 0{r^). (3-26) 

2i53 

We find that both gravitational potential and gravitational force approach constant 
values, which is consistent with the observation in Fourier space that effective grav- 
itational constant decays like oc A;^^ at large k. 

It is important to notice that in the limit Ai — )• 1, ^2 — 0, Si — )• 1, -B2 — ^ and 
i?3 — 7- one recovers the known GR result, that '\s Q = 1. Although we want to be 
as much general as possible in this study, we notice that this happens for theories 
for which Fn^Fg^^g all go to at early times. On the other hand, all these terms 
are supposed to be of the same order at late times. At the typical scales at which 
linear perturbation theory can be used to test a theory, A;/(ao-f^o) ~ 100, one can set 
some constraints as follows 

1. I^s/^zl X 10^ < 1. This case happens when the modifications of gravity are 
either of the f{R) kind or very weak even today, or when the expansion of 
the universe is almost de Sitter {H ~ 0). In particular the theory reduces 
to the special cases (for which Fg — Fji^g 0) discussed later on. This case 



*•* This ansatz does not solve Eq. (|3-28[1 . However, we can give it as an initial condition for 
together with Sm ~ 0. The fact that for large k Q ex implies that Sm will not evolve rapidly but 
only in Hubble time. In other words, we can consider it as a good approximation in a short time 
range. 
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although possible is not interesting enough, because one would need some other 
experiment to disentangle these kind of theories from, say, the f{R) ones. 
2. I-B3/S2I X 10'' < 1 and \B2/Bi\ x 10'' < 1. When this case happens, since A2 
is typically the same order of magnitude as B2 unless the model parameters 
are fine-tuned, in general we also have IA2M1I x 10'' < 1. This case leads to a 
matter spectrum close to GR. 
It is then evident that when any of the previous two conditions are not satisfied we 
are in the full Modified Gravity regime, for which bounds can be set. In particular 
for scales at which dominates, since Q — )• we would expect the matter power 
spectrum to reduce, as gravity is weaker at those scales. On the other hand, at 
scales for which B2 is the dominant term, the quantity one needs to extract from the 
particular model at hand \s Q ^ A2/ B2. If this number is larger or smaller than 1 /F 
then one obtains an enhancement or a suppression of the matter power-spectrum. 
Finally, in general one may find both these two different behaviors, which depend 
on the scale. 

3.4. Evolution equation for the matter density perturbation 

From (|3-12|) and (j3-13p . the following evolution equation for 5m is derived: 

5m + 2H5m + ^^ = 0. (3-27) 



Putting the Poisson equation (|3-7p into this equation yields a closed evolution equa- 
tion for 5m, 

5m + 2H5m - AnGNQpmSm = 0. (3-28) 

3.5. Parameter E 

In this subsection, we provide expressions of a parameter S that are sensitive 
to weak lensing for f{R, Q) theories. This parameter is defined by 

r^la+ry), (3-29) 

where rj is defined by 

7? = -|. (3-30) 

The photon propagation is sensitive only to a so-called lensing potential <P — ^ 
and the parameter S appears when we rewrite <P — ^ m. terms of 5m- Therefore, E 
parametrizes the relation between the lensing potential and the density perturbation. 
We find that E is given by 



Ci + C2{k/aY 
Bi+B2{k/af + B?Xk/ay 



^= ^ . :\:r.r;\,. (3-31) 



where Ci and C2 are given by 



Ci = F, 



(3-32) 
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Co 



Fr + 4Fcj(2/72 + H) + 16H^^g{H^ + H) 



+ 8H {Fg + A{H^ + H)Cg 



(3-33) 



Exactly as with the effective gravitational constant, U approaches a constant at large 
distance and decays in proportion to at short distance. 

§4. Special cases 

The consequences for generic f{R,G) theories derived in the last section do not 
apply if the model satisfies a condition, 



{Fnig - Fl) = 0. 



(4-1) 



For example, f{R) and R + f{Q) theories belong to this class. However, there are 
infinite number of other theories which satisfy this condition.-^ If this condition is 
satisfied, vanishes exactly. 

Interestingly, ()4-ip is exactly the same condition as the one for the absence of fc^- 
term in the dispersion relation for the fast oscillating modes.^ If the theory satisfies 
()4-ip . the dispersion relation becomes ojI = Ak'^, where explicit form of A depends on 
the theory. For example, for f{R) theories, A = 1/a?. Therefore, the fast oscillating 
modes propagate with a velocity of light and the sub-horizon approximation we have 
used can be still applied. For R + f{Q) theories, however, it was found inP that 
modes corresponding to the fast oscillating ones become highly unstable on small 
scales in the radiation/matter dominated era. In such a case, the unstable modes 
cannot be neglected and hence we should include them in the perturbation analysis. 
In this paper, we have assumed that modes corresponding to the fast oscillating ones 
in the theory satisfying ()4Tp are stable and safely decouple from the other modes. 

For the special case where (j4-ip holds, the effective gravitational constant ap- 
proach finite values for large k. Decay of Q that we have observed in the last section 
does not occur for the special case. Correspondingly, the gravitational potential 
p-24p around a mass point in this case is given by 



GnM Ai 



1 + 



B1A2 - A^ B 



1-D2 



AIB2 




■exp -\ — 



(4-2) 



The second exponential term acts as shifting the gravitational constant by a fac- 
tor B1A2/ {A1B2) as we go down to the smaller scales. This exponential type of 
gravitational potential was also derived for f{R) theories.'^"^ 

The power-law decay at short distance is also absent for E. We see that it 
asymptotically approaches C2/B2 at short distances. In the following, just for the 
purpose of demonstration, we will give expressions of Q and S for f{R) and f{R + 
Q/A^) theories. 
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4.1. Special case I : f{R) theories 

Let us first consider f{R) theories. In this case, we have 

Q = li±l^, (4-3) 

17=1. (4.4) 

Therefore, we have U = Q at large scales and U = at short scales. 

4.2. Special case II : f{R + Q j theories 

Let us next consider f{R + Q j theories, where yl is a constant of mass di- 
mension. We can easily check that this model satisfies ()4-ip . In the limit /l — )• oo, 
the theory reduces to the f{R) theory. 

In this model, and [i are given by 



F yl6i7+(yl2 + 4^2)|y^|2 



_ yl2 yl4 + (yl2 + 4^^ + 4ij)(3yl2 + 127/2 + g^):^^ 

yl6F+(yl2+4i/2)VF^g ' ^^'^^ 

where PV^ is defined by 

W = mHF{A^ + 4i?2 + 4ij) + A^F{M^ + 12i?2 + I6ij) - 4(yl2 + 4^^)^. (4.7) 

Therefore, we have i7 = Q at large scales and 

^ ^ 3^1^+12^2 + 8^ 

4(yl2 + 4/72 + 4ij)^' ^ ^ 

at short scales. 



§5. Discussion and Conclusion 



We have shown that in f{R,Q) theories, the effective gravitational constant 
generically decays in proportion to k~'^ on sub-horizon scales. This means that the 
gravitational force between the two massive objects is weaker than the one in GR 
and such a deviation will be larger on smaller scales. Then, we naively expect that 
the strongest constraint on f{R,Q) theories will be obtained either from the solar 
system constraints or Cavendish type experiments on the Earth. However, the direct 
application of our result to such small scale systems is not justified since the back- 
ground spacetime is no longer FLRW universe. Instead, the background spacetime 
we should use for the solar system, for example, is the spherical symmetric static 
metric that must be close to the Schwarzshild metric. Just to illustrate this issue 
more in detail, we will consider here, as one example, the solar-system constraints 
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on the dark energy modelP^ 

f {R, g) = R + X-^arctan ( ^] - aX^/g'^, (5-1) 

where a, A and are positive constants, was studied (many other models were 
also analyzed inP^). The Gauss-Bonnet correction to GR is highly-suppressed in 
the solar-system by e = \fQll\fQ's ^ 1 where Qa is the Gauss-Bonnet term at the 
solar-system. It was shown that the spherical-symmetric metric: 

= -^(r)dt2 + -l-dr2 + r^dl?^ (5-2) 
B[r) 

can be expanded in e as, 

A{r) = 1 - !£ + eci f-V + 0{^), B(r) = 1 - !l + f-V + 0(e')(5-3) 

where ci,C2,p and q are model-dependent 0{1) constants. 

For this kind of background, it is clear that the effective gravitational constant 
is close to the Newton's constant and we do not see the /c~^ scaling. Therefore, 
the diminishing of the gravitational constant we have observed occurs only when 
the objects whose background spacetime can be treated as FLRW universe, that 
is, when the expansion ()5-3p breaks down 0. For some models, the distance r at 
which the expansion breaks down can be much smaller than the horizon scale (But 
it must be still larger than the solar-system size for such a model to be viable.). 
In such a case, there may be a region of distance scales where the weakening of 
gravity may be significant and its signature may be imprinted in the patterns of the 
cosmic structures such as group of galaxies as less structures than GR, which can be 
tested by observations. In principle, we can construct f{R,Q) theories that possess 
those properties by requiring that /(i?, G) becomes very close to GR for large R and 
deviates from GR when R and Q are 0{Hq) and 0{Hq), where Hq is the Hubble 
constant. 

We also found that if the theory satisfies H^S = 0, then the weakening of gravity 
does not occur. Interestingly, exactly the same combination appears in the dispersion 
relation for the fast oscillating modes. We conclude that the short scale behavior 
of the perturbations and the gravity are qualitatively dependent on whether H'^S 
vanishes or not. If, at the scales where perturbation theory can be applied to study 
the data, this quantity does give some non negligible contribution, then we expect 
the matter spectrum, at the same scales, to be in general suppressed. 
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